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{2L|' Abstract 

. We present numerical results for the two loop QCD corrections to the ghost gluon 

vertex, the quark gluon vertex, and the triple gluon vertex in the massless limit at 
the symmetric point. We restrict ourselves to the tensor structures existing in the 
tree level. The corrections are used to examine different coupling constant definitions 



5^ I in momentum subtraction schemes and the corresponding three-loop /? functions. 



1 Introduction 

There is no unique definition of the QCD couphng constant - a^: its value depends on the 
renormahzation prescription employed. Within pQCD the definition which is most often 
used is based on the MS-scheme [1, 2]. Such a definition is of great convenience for dealing 
with inclusive physical observables dominated by short distances (for a review see [3]). 
On one side the underlying use of dimensional regularization makes advanced calculations 
possible: e.g. the corresponding /J'^^-function is available now with four-loop accuracy [4]. 
On the other one, this makes the MS-scheme virtually inapplicable in cases when one 
chooses to utilize a different regularization, including, most regrettably, the lattice one. In 
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addition, the physical meaning of the MS normahzation parameter // is not transparent and 
leads to the well-known ambiguities when considering the decoupling of heavy particles. 

These shortcomings are absent for a wide class of so-called momentum subtraction 
(MOM) schemes. The MOM-schemes require the values of properly chosen Green functions 
with predefined /^-dependent configurations of external momenta to be fixed (usually to 
their tree values) independently on the considered order. Practical calculations can then 
be performed with any regularization, including the lattice one. 

Unfortunately, even for massless QCD, there are infinitely many possibilities to define 
a momentum subtraction renormalization scheme. Not only is there an ambiguity which 
vertex to subtract, and at which exact configuration of external momenta. In addition, 
there is the freedom to use a certain linear combination of the scalar functions appearing 
in the gluon and quark vertices, which can be related to fixed polarization states of the 
external particles [5]. 

Recently, the momentum subtraction approach has been heavily used to relate lattice 
results for quark masses [6, 7, 8] and coupling constants [9, 10, 11, 12, 13, 14, 15] to their 
perturbatively determined MS counterparts. 

In one of these papers [13] it has been argued that the knowledge of three-loop coeffi- 
cients for the corresponding /5-functions is necessary and oven the four-loop contributions 
should be taken in to account. This is because the accessible energy ranges in these cal- 
culations are just reaching a level where perturbative QCD calculations start to be valid 
approximations . 

The results of these lattice calculations have mainly been analyzed within the so- 
called MOM-scheme. The scheme is defined using an asymmetric subtraction at qf = 
?2 — ~A*^)?3 = 0- The perturbative calculations with such a definition are certainly 
much simpler than those for the schemes employing the symmetric subtraction point 
Qi — Q2 — Qs — ~A*^- ^ consequence only for that scheme the non-trivial three-loop 
coefficient of the /9-function is known^. On the other hand, on general grounds it is rather 
clear that defined through the symmetric subtraction point should be less prone to all 
kinds of non-perturbative efi:ects than the one based on the asymmetric subtraction. 

In this paper we suggest a new (approximate) way to compute (massless) vertices at the 
symmetric point by using a large momentum expansion. The method reduces the problem 
to evaluation of massless propagators. We compute the two loop QCD corrections to the 
ghost gluon vertex, the quark gluon vertex, and the triple gluon vertex in the massless 
limit at the symmetric point in the MS-scheme. The obtained results are then used to 
construct different coupling constant definitions in momentum subtraction schemes and 
the corresponding three-loop P functions. 

This letter is organized as follows: In the next Section we discuss the crucial difference 
concerning possible non-perturbative corrections to the running of the coupling constant 
defined with subtraction at the symmetric momentum configuration (which is used in 
MOM-schemes) and the asymmetric one (used in the so-called MOM-schemes). Section 3 
briefiy describes our calculational method. In Section 4 the results for the vertex corrections 
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at two loops are given. Then we derive the relations between the coupling constants in the 
MOM- and the MS-scheme (Section 5) and close by a short discussion on the effects of our 
results on the three loop j3 function in the MOM-scheme (Section 6). 

2 MOM vs. MOM 

A well-known way to take into account at least some of non-perturbative physics is the 
QCD sum rules method [17] based on the use of Operator Product Expansion (OPE) (for 
a recent review, see [18]). For instance, let us consider a correlator of two local operators 

G^'^iq) = I j e'''"da;(0|T{A(a;)5(0)}|0). (1) 

At large (and Euclidean) external momentum transfer q the correlator (1) can be 
schematically represented in the form 

^^^(^) 1^ +ECn^^(?)(On(0)), (2) 

n 

where the first term, G'^'^(g), stands for the purely perturbation theory contribution while 
the sum describes the factorization of the non-perturbative effects due to the large distances 
(hidden inside of vacuum expectation values (VEV) of various composite operators) and 
the coefficient functions C^^{q). The latter correspond to the short distance (of order 
l/|g'|) contributions coming from the integration with respect to x in (1) and, thus, are 
computable within perturbation theory. 

Next, let us consider a three-point correlator 

G^^^{p, Q)^i^ J e''^''+'Pydxdy{0\T{A{x)B{0)C{y)}\0) (3) 

in the kinematical regime relevant for the MOM-scheme, that is for both q and p being 
large. In this case one can straitforwardly apply the OPE for three operators and write: 

G^^'^'iP^^i) Gr^(p,g) + ECn'^''(P,?)(On(0)). (4) 

Here the coefficient functions C^^^{p, q) are computable as an expansion in the (small) 
coupling constant as{[j) normalized at (large) momentum scale of order \p\ ^ \q\. An 
important feature of both Eqs. (2) and (4) is that the VEV of composite operators are 
universal and do not depend on the correlator. 

The situation is significantly different for the case when only one external momentum, 
say, q is large and another one is fixed or even set to zero (the latter case corresponds to 
a MOM-scheme). Indeed, in this case the representation (4) should be replaced with [19] 

i 

+ E^n''(5) ^ / e^^M|/(0|T{C(|/)O.(0)}|0). (5) 
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In this representation the first sum comes from the integration region where both x and y 
are small and of order l/\q\. Note that the functions C^^'^{p,q) are in general different 
from C^^'^{p, q). 

The second sum results from integration regions of small x ~ '^/\q\ only (where the OPE 
of two operators A and B can be employed). As the momentum p is assumed to be small or 
even zero the two-point correlators in the second sum can not be perturbatively computed 
and should be considered as phenomenological quantities analogous to " condensates" . 

Thus, wc conclude that the three-point correlators employed in MOM-schemes contain, 
in in addition to the VEV of composite operators, an extra source of the non-perturbative 
corrections — VEV of bifocal operators — in comparison to the same correlators with 
symmetrical pattern of the external momenta. 

3 Method of the calculation 

The analytic computation of arbitrary two loop three point functions seems to be exchidcd 
at the time being. Therefore our computation is based on the method of asymptotic 
expansions [20] which reduces the complexity of the integrals. We consider three point 
graphs with incoming momenta gi, ^2, and — (q'l + ^2) and identify one momentum as large, 
say qi. The ordinary large momentum procedure is used to find an expansion with respect 
to and qi.q2/qi- The symmetric point further demands that {qi + ^2)^ = qf, so qi.q2 

can be reexpresscd by Setting = zqf yields a series in the variable z. 

Using the described method we found very good numerical agreement with the analyt- 
ical results for the one loop vertex corrections given in [21, 22] and the one and two loop 
scalar three point integrals given in [23] by computing only 4 to 6 terms of the expansion. 

The use of the large momentum procedure results in an asymmetric approach and 
marks one momentum. Therefore one obtains at least two independent series by choosing 
different momentum distributions. The numerical values of the series must coincide at the 
symmetrical point {z = 1). This provides a powerful check on our method and can be 
tested by examining the non-symmetric scalar ladder type topology computed in [23]. 

Wc use the scries with the best convergence to estimate the mean value and the variation 
of independent expansions as an indication for the size of the uncertainty. This description 
is in most cases equivalent to an error estimation based on the size of the highest expansion 
term available. 

Throughout the whole computation we made use of several computer programs. The 
diagrams were generated using QGRAF [24] and the large momentum procedure was applied 
by EXP [25] . The actual evaluation of the integrals was done using the MINCER [26] package 
written in FORM [27] . 

4 Two loop vertex corrections 

We parametrize the QCD vertices at the symmetrical point {pf = pi = pi = p^) as follows: 
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f«(-/.2) 




32 ^ 96 




AW(-//2) 


2 ~ 12 


13 13 T 
16 96 




rW(V) 




3 _|_ 23 7- 
32 288 
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Table 1: One loop results — by SU(N) color factors. 







CfCa 


CI 


CATuf 


CpTuf 








0.197(1) 


-0.151(2) 




A(2)(-;i2) 


0.206(4) 


-0.20(4) 


0.679(1) 


-0.4968(4) 


-0.0211(4) 


r(2)(-/.2) 






-0.22(4) 


0.65(7) 


-0.408(10) 



Table 2: Two loop results — by SU(N) color factors. 



• ghost gluon vertex (pi is the incoming momentum of the out-ghost) 

^t''{pi,P2,P3) = igsr'y; {g,.r{p') + •••), (6) 

• quark gluon vertex 

Al,,{puP2,P3) = 9sTt, (7mA(/) + ■■■), (7) 

• triple gluon vertex 

-la&c / 



KMpi^p^^p^) = -^9sr"'x 



{gijiv {pi - P2)x + 9v\ ip2 - P3)^ + g\,x iP3 - Pi)^) r(p^) h — ] (8) 

with the strong coupling constant Qg, the structure constants of the SU(N) Lie algebra 
jabc^ and the matrices of the fundamental representation T^j. The terms contained in the 
ellipsis are defined in such a way that they do not interfere with the displayed form factor. 
This means that e.g. for Eq. (7) the dots could contain only the structures Pi,i/iPj and 
so on. Note, that our choice of the kinematical structures (6-8) is, probably, the simplest 
one. However, our approach could be applied, if necessary, to any other choice, should the 
latter be more convenient and/or natural for a specific physical problem. 

With Eqs. (7) and (8) we follow the prescription of [21]. Concerning the ghost gluon 
vertex our definition partly avoids the ambiguity connected with the identification of the 
vertex form factors and differs from the one used in [21]. 
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Each form factor receives quantum corrections: 

r = i + ^rw + f^yr(^) + ... (9) 

TT V TT / 

and similarly for T and A. 

In the following we restrict ourselves to Landau gauge so the gluon propagator is taken 
to be purely transversal. A further simplification can be achieved by observing that the 
MOM-scheme is usually defined at = —fi^ with the 't Hooft mass fi. Therefore all 
logarithms log(— drop out. All quantities displayed in the following are computed 
using the MS-scheme for the renormalization. 

The one loop corrections to the ghost gluon vertex, the quark gluon vertex, and the 
triple gluon vertex can be found in [21]. There the integral 

7 = -2 ^^^^^^^ = 2.3439072 . . . (10) 
Jo — X + 1 

was defined. The results are displayed in table 1 where we have distinguished between 
contributions proportional to different color factors of the SU(N) (for SU(3) Cp — 4/3, 
Ca = 3, and T = 1/2, Uf represents the number of fermions). 

Our numerical results for the two loop contributions are summarized in table 2. The 
number in the brackets denotes our error estimation on the last digit. 

Inserting the values of the color factors for SU(3) and hnearly adding the errors one 
would surely overestimate the uncertainty. Therefore we first computed the series with 
color factors already inserted and then extracted the results which read 

f(2)(-y^2^ = 1.770(9) - 0.1727(3) n/, 
A('^\-I_i^) = 5.590(6) - 0.7594(2) n/, 

r^^\-l_i^) = -2.0(4) + 0.72(9) n/. (11) 



5 Relations between coupling constants 

The MOM renormalization condition implies that the value of one of the above defined 
functions should be equal to 1 when the corresponding MOM-scheme is used. This allows 
for the computation of the coupling constant in that MOM-scheme (a^^^^) as a perturba- 



tive series in the coupling constant of the MS-scheme 



a 



MS^ 



which is of the form 



a 



MOM 



a 



MS 



1 + ^ 
TT 



— ^ 1 (d2o + d2inf + d22nfj + ■ 



;i2) 



We will refer to MOM quantites computed on the basis of the triple gluon vertex, the quark 
gluon vertex, and the ghost gluon vertex as MOMggg, MOMq, and MOMh, respectively. 

Using the well known results for the self energies of the quark, the ghion, and the ghost 
at one and two loops in the MS-scheme we find the coefficients shown in table 3. 
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dio 


du 


d20 


1^21 


C^22 


MOMh 


49 _|_ 1 r 
12 16 


5 

18 


35.88(2) 


-5.0707(6) 


25 
324 


MOMq 


89 85 T 
16 144-' 


5 
18 


29.53(1) 


-5.1961(4) 


25 
324 


MOMggg 


11 I 23 7- 
2 48 


_1 _ 2j 

3 9 


59.8(8) 


-12.6(2) 


432 ~ 54-' ~ 81 



Table 3: Coefficients of tfie relation between the coupling constants in the MOM- and the 
MS-scheme. 

6 Beta function 

We define the perturbative series of a generic (3 function by 

,.=A— = /'(«.) = -f-)'Eftf-)'- (13) 



i>0 



TT 



The P function in the MOM-schcmc (Z?^^*-"^) can be computed directly from the results 
presented in the last section using the relation 



J MOM/ 2\ dn^^^ 

^MOMj-^MOMj _ ^2 ^ ) _ pMS^^s 



dA(2 TT " da^s 



(14) 



In our framework the first two terms in the perturbative series of the /? function are 
scheme-independent. To compute Z?^*-*^^ the knowledge of the two loop vertex corrections 
is sufficient if is given. 

^MS jg j^nQ^n fQ^p loops [28, 29, 30, 4] and the first three terms read: 

,r ^ i(T--H- 

1 /34_o „ 20 



= ^{^C% + 2Cf.Tnf-^CFCATnf-^CiTnf 

+ fc,T^n}+^-^CAT^n}). (15) 

The vertex corrections given above lead to the following values for the three loop /3 
function in the MOM-scheme: 



1 /2857„, „„o„ 205 1415 



^MOMh _ 44.82(5)- 9.730(5) n/ + 0.3276(1) nl 



n 



/' 



^MUiviq ^ 28.86(3)- 9.206(3) + 0.35322(7) 
^MOMggg ^ 24(2) + 0.04(63)71/ - 1.05(3) nj + 0.0415330 nj. (16) 
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Note, that a large uncertainty of the ^/-coefficient for the MOMggg-scheme is the result 
of an accidental cancellation between (3^^ and the corrections induced by the one and two 
loop vertex corrections in this case. 

It is instructive to compare the MOMggg-scheme with the (standard) MS-scheme as well 
as with the MOMgg one. The latter is defined with the help of the asymmetric subtraction 
of the triple gluon vertex (more details in [11, 16]). The three- loop contributions to the 
corresponding /3-functions read 

/5p = 22.3203 -4.36892n/ + 0.0940394 nj, 
^MOMgg ^ 37.6899 - 5.57013n/- 0.223177n5 + 0.0138889 nj. (17) 

Thus, one observes that (at three loops) the MOMggg-scheme is numerically signifi- 
cantly closer to the MS-scheme that to the MOMgg one. 

Very recently, in [31] an attempt has been made to "predict" the value of (32^'~'^^\nf=o 
with the help of an OPE analysis of the flavourless non-perturbative gluon propagator and 
the symmetric triple gluon vertex in the Landau gauge. They have found 

^MOMgggi r/q\ /oMOMggi 

while our result is 

^MOMggg|^^^^ = 0.64(5) °^^^U,=o. (19) 

We believe that the rough (that is in the sign and overall magnitude) agreement between 
Eqs. (18) and (19) could be improved by taking into account anomalous dimensions of 
power corrections as discussed in [31]. 



Conclusion 

We have presented a new approach to approximately compute three-point Green functions 
in massless theories. The approach is based on the use of asymptotic expansion methods. 
It reduces the problem of computing of a three-point N-loop correlator (with Euclidean 
external momenta) to the calculation of the N-loop massless propagators. We have checked 
the method by recomputing a few analytically known results. In all cases the resulting series 
proved to be accurate approximations to the exact results. 

We have applied the approach to compute the two loop QCD corrections to the ghost 
gluon vertex, the quark gluon vertex, and the triple gluon vertex in the massless limit 
at the symmetric momentum point. The results have been used to construct the corre- 
sponding three-loop (3 functions. In principle, the calculation could be extended to include 
one more loop (as the three- loop massless propagators are still accessible with MINCER). 
Unfortunately, for the time being the hardware constraints look insurmountable. 



8 



Acknowledgements 

We would like to thank J.H. Kiilin for numerous fruitful discussions and A.I. Davydychev 
for providing us with information about the results of [22] . We are grateful to the authors 
of the work [31] for sending us a draft version of it. 

This work was supported by DFG under contract FOR 264/2-1, the Graduiertenkol- 
leg "Elementarteilchenphysik an Beschleunigern" and the DFG-Forschergruppe "Quanten- 
feldtheorie, Computeralgehra und Monte-Carlo-Simulationen" . 

References 

[1] G. 't Hooft, Nucl. Phys. B 61 (1973) 455. 

[2] W.A. Bardeen, A.J. Buras, D.W. Duke, and T. Muta, Phys. Rev. D 18 (1978) 3998. 
[3] K.G. Chetyrkin, J.H. Kiihn, and A. Kwiatkowski, Phys. Reports 277 (1996) 189. 
[4] S.A. Larin, T. van Ritbergen, and J.A.M. Vermaseren, Phys. Lett. B 400 (1997) 379. 
[5] E. Braaten and J.R Leveille, Phys. Rev. D 24 (1981) 1369. 

[6] G. Martinelh, C. Pittori, C.T. Sachrajda, M. Testa, and A. Vladikas, Nucl. Phys. B 
445 (1995) 81. 

[7] E. Franco and V. Lubicz, Nucl. Phys. B 531 (1998) 641. 

[8] D. Becirevic, V. Gimenez, V. Lubicz, and G. Martinelli, Phys. Rev. D 61 (2000) 
114507. 

[9] D. Henty, C. Parrinello, and C. Pittori, Proceedings of the International Europhysics 
Conference on High Energy Physics, Brussels (1995), 239. 

[10] B. AUes, D. Henty, H. Panagopoulos, C. Parrinello, C. Pittori, and D. Richards, 
Nucl. Phys. B 502 (1997) 325. 

[11] P. Boucaud, J. P. Leroy, J. Micheh, O. Rene, and C. Roiesnel, JHEP 12 (1998) 004. 

[12] P. Boucaud, J. P. Leroy, J. Micheh, O. Rene, and C. Roiesnel, JHEP 10 (1998) 017. 

[13] D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Rene, J. Rodriguez-Quintero, and 
C. Roiesnel, Phys. Rev. D 60 (1999) 094509. 

[14] D. Becirevic, P. Boucaud, J.R. Leroy, J. Micheli, O. Rene, J. Rodriguez-Quintero, and 
C. Roiesnel, Phys. Rev. D 61 (2000) 114508. 

[15] J. I. SkuUerud, UKQCD col, Nucl. Phys. Proc. Suppl. 63 (1998) 242. 



9 



[16] K.G. Chetyrkin and A. Retey, preprint hep-ph/0007088. 

[17] M.A. Shifman, A.I. Vainstein, and V.I. Zakharov, Nud. Phys. B 147 (1979) 
385,448,519. 

[18] F. V. Gubarev, M.I. Polikarpov and V.I. Zakharov, Nud. Phys. Proc. Suppl. 86 (2000) 
437. 

[19] I.I. Balitsky and A.V. Yung, Phys. Lett. B 114 (1982) 53; 

B.L. loffe and A.V. Smilga, Pis'ma Zh. Exp. Teor. Fiz. 37 (1983) 250; 

K.G. Chetyrkin, S.G. Gorishny, A.B. Krasulin, S.A. Larin, and V.A. Matveev, preprint 

INRP-0337 (1984). 

[20] see e.g., V.A. Smirnov, Mod. Phys. Lett. A 10 (1995) 1485, and references therein. 

[21] W. Celmaster and R.J. Gonsalves, Phys. Rev. D 20 (1979) 1420. 

[22] A.I. Davydychev, R Osland, and O.V. Tarasov, Phys. Rev. D 54 (1996) 4087, erratum 
ibid. D 59 (1999) 109901. 

[23] N.I. Ussyukina and A.I. Davydychev, Phys. Lett. B 332 (1994) 159. 

[24] P. Nogueira, J. Comp. Phys. 105 (1993) 279. 

[25] T. Seidensticker, Diploma thesis (University of Karlsruhe, 1998), unpubhshed. 

[26] S.A. Larin, F.V. Tkachov, and J.A.M. Vermaseren, preprint NIKHEF-H/91-18 (1991). 

[27] J.A.M. Vermaseren, Symbolic Manipulation with FORM, CAN (1991). 

[28] D.J. Gross, F. Wilczek, Phys. Rev. Lett. 30 (1973) 1343; 
H.D. Politzer, Phys. Rev. Lett. 30 (1973) 1346; 

G. 't Hooft, report at the Marseille Conference on Yang-Mills Fields, 1972. 

[29] W.E. Caswell, Phys. Rev. Lett. 33 (1974) 244; 

D. R.T. Jones, Nud. Phys. B 75 (1974) 531; 

E. S. Egorian, O.V. Tarasov, Teor. Mat. Fiz. 41 (1979) 26. 

[30] O.V. Tarasov, A.A. Vladimirov, and A.Yu. Zharkov, Phys. Lett. B 93 (1980) 429; 
S.A. Larin, J.A.M. Vermaseren, Phys. Lett. B 303 (1993) 334. 

[31] P. Boucaud, A.L. Yaouanc, J. P. Leroy, J. Micheh, O. Pene, and J. Rodriguez-Quintero, 
preprint hep-ph/0008043. 



10 



